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Abstract — [] In this paper, achievable rates of imperfectly- 
known fading relay channels are studied. It is assumed that 
communication starts with the network training phase in which 
the receivers estimate the fading coefficients of their respective 
channels. In the data transmission phase, amplify-and-forward 
and decode-and-forward relaying schemes are considered, and 
the corresponding achievable rate expressions are obtained. The 
achievable rate expressions are then employed to identify the 
optimal resource allocation strategies. 

I. Introduction 

In wireless communications, deterioration in performance is 
experienced due to various impediments such as interference, 
fluctuations in power due to reflections and attenuation, and 
randomly-varying channel conditions caused by mobility and 
changing environment. Recently, cooperative wireless com- 
munications has attracted much interest as a technique that 
can mitigate these degradations and have the performance 
approach to the levels promised by multiple-antenna systems. 
Cooperative relay transmission techniques have been studied 
in [1] and [2] where several two-user cooperative proto- 
cols have been proposed, with amplify-and-forward (AF) and 
decode-and-forward (DF) being the two basic modes. In [3], 
three different time-division AF and DF cooperative protocols 
with different the degrees of broadcasting and receive collision 
are studied. In general, the area has seen an explosive growth 
in the number of studies (see e.g., [4], [5], [6], [7], [8] and 
references therein). However, most work has assumed that 
the channel conditions are perfectly known at the receiver 
and/or transmitter sides. Especially in mobile applications, this 
assumption is unwarranted as the randomly-varying channel 
conditions can be learned by the receivers only imperfectly. 
Recently, Wang et al. in [9] considered pilot-assisted trans- 
mission over wireless sensory relay networks, and analyzed 
scaling laws achieved by the amplify-and-forward scheme in 
the asymptotic regimes of large nodes, large block length, and 
small SNR values. In this study, the channel conditions are 
being learned only by the relay nodes. 

In this paper, we study the achievable rates of imperfectly- 
known fading relay channels. A priori unknown fading coeffi- 
cients are estimated at the receivers with the assistance of pilot 
symbols. Following the training phase, AF and DF relaying 
techniques are employed in the data transmission. Achievable 
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rates for these schemes are used to find the optimal resource 
allocation strategies. 

II. Channel Model 

We consider the three-node relay network which consists 
of a source, destination, and a relay node. Source-destination, 
source-relay, and relay-destination channels are modeled as 
Rayleigh block-fading channels with fading coefficients de- 
noted by h sr , h sd , and h rd , respectively, for each channel. Due 
to the block-fading assumption, the fading coefficients h sr ~ 
CAA(0, a sr 2 ), h sd ~ CN(0,a sd 2 ), and h rd ~ CA/"(0, u rd 2 ) El 
stay constant for a block of m symbols before they assume in- 
dependent realizations for the following block. In this system, 
the source node tries to send information to the destination 
node with the help of an intermediate relay node over the 
coherence block of m symbols. The transmission is conducted 
in two phases: network training phase and data transmission 
phase. Over these phases the source and relay are subject to 
the following power constraints: 

||x M || 2 + £{||x s || 2 }<mP s , (1) 

Hx^H 2 + £{||x r || 2 } <mP r . (2) 

where x s t and x r jt are the source and relay training signal 
vectors respectively, and x s and x r are the corresponding data 
transmission vectors. 

A. Network Training Phase 

Each block transmission starts with the training phase. In the 
first symbol period, source transmits a pilot symbol to enable 
the relay and destination to estimate channel coefficients h sr 
and h s d- In the average power limited case, sending a single 
pilot is optimal because instead of increasing the number of 
pilot symbols, a single pilot with higher power can be used. 
The signals received by the relay and destination, respectively, 
are 

y r ,t = h sr x St t + n r , (3) 
Vd,t = h sd x a<t + n d . (4) 

2 x ~ CAf(d, a 2 ) is used to denote a proper complex Gaussian random 
variable with mean d and variance a 2 . 



Similarly, in the second symbol period, relay transmits a 
pilot symbol to enable the destination to estimate the channel 
coefficient h rd . The signal received by the destination is 



Vd t = KdXrd + n d . 



(5) 



In the above formulations, n r ~ C7V(0, N ) and n d ~ 
CAf(0, No) represent independent Gaussian noise samples at 
the relay and the destination nodes. 

In the training process, it is assumed that the receivers 
employ minimum mean-square error (MMSE) estimation. Let 
us assume that the source allocates S s of its total power for 
training while the relay allocates 5 r of its total power for 
training. As described in [12], the MMSE estimate of h sr is 
given by 



a; r \/S s mP s 
a 2 r 8 s mP s + N 



Vr,t, 



(6) 



where j/ nt ~ £/V(0, a 2 r d s mP s + N ). We denote by h sr the 
estimate error which is a zero-mean complex Gaussian random 
variable with variance 

var(li sr ) = 9 r~F, , », ■ (') 



a 2 sr 5 s mP s + A^o ' 



Similarly, we have 



- a 2 sd ^J5 s mP s 
a 2 d 5 s mP s + N 
y dtt ~CM(0,o- 2 sd 6 s mP s + No), 



var(h sd ) 



a 2 d 5 s mP s + N Q ' 



(8) 



(9) 



(J 2 d \J8 r mP r 



<J 2 rd 5 r mP r + No****' 
y r dtt ~ CAf(0,a 2 rd 6 r mP r + No), 

rl \ a rd N o 
var(h rd ) - 



a; d S r mP r + No 



(10) 
(ID 



With these estimates, the fading coefficients can now be 
expressed as 

(12) 



hrd h r d h r d- 



(13) 
(14) 



B. Data Transmission Phase 



The practical relay node usually cannot transmit and receive 
data simultaneously. Thus, we assume that the relay works 
under half-duplex constraint. We further assume that the relay 
operates in time division duplex mode. As discussed in the 
previous section, within a block of m symbols, the first two 
symbols are allocated for channel training. In the remaining 
duration of m — 2 symbols, data transmission takes place. 
First, the source transmits an (m — 2)/2-dimensional symbol 



vector x s which is received at the the relay and the destination, 
respectively, as 

y r = /j si .x s + n r , (15) 



y d = h sd x s + n d . 



(16) 



Next, the source becomes silent, and the relay transmits an 
(to — 2)/2-dimensional symbol vector x. r which is generated 
from the previously received y r [1] [2]. This approach cor- 
responds to protocol 2 in [3], which realizes the maximum 
degrees of broadcasting and exhibits no receive collision. 
Thus, the destination receives 

y r d = h rd x r + n d . (17) 

After substituting ([T2j-([T4]> into (Q3), (O, ([17]) we have 

(18) 

(19) 

(20) 



y d = h sd x s + h sd x s + n d , 
y r d = h. rd yi r + hrdXr + n d . 



The input vectors x s and x r are assumed to be composed of 
independent random variables with equal energy. Hence the 
corresponding covariance matrices are 

to — 2 

+ 1 2(1 -<UmP r 
£{x r xt} = -i V^I, (22) 

to — 2 

where I is the (to — 2)/2-dimensional identity matrix. 

III. A Capacity Lower-Bound For AF 

In this section, we consider the AF relaying scheme and 
calculate a capacity lower bound using similar methods as 
those described in [11]. The capacity of the AF relay channel 
is the maximum mutual information between the transmitted 
signal x s and received signals y d and y d given h sr , h sd , h rd . 
Thus, the capacity is 

C = sup — 7(x s ;y d ,y^|/i sr ,^ S£; ,/i rd ). (23) 

Note that this formulation presupposes that the destination has 
the knowledge of h sr . Hence, we assume that the value of h sr 
is forwarded reliably from the relay to the destination over 
low-rate control links. 

Our method for finding a lower bound obtains h sr , h sd and 
h rd , relegates the estimation error of channel estimates to the 
additive noise, and then considers only the correlation (and 
not the full statistical dependence) between the resulting noise 
and the transmitted signal. We then obtain a lower bound by 
replacing the resulting noise by the worst case noise with the 
same correlation. Let us assume that 



z d = h sd x s + n d , 
z d = h rd x r + n d , 



(24) 
(25) 
(26) 



2 



are the noise vectors which has the following covariance 
matrices: 



£{ ZrfZ t} = <7? st £{x s xt} + iVoI, 
E{z r d z r J} = o? E{xr4} + N I. 



(27) 
(28) 
(29) 



where i = 3,4, (m — 2)/2 + 2. With the above notation, 
we can write the input-output mutual information as 

(m-2)/2+2 

I(-x. s ]y d ,y r d \h sr ,'h sd ,h rd ) = V /(x s [i];yd[i]|/i. SI .,^ s d,/i r rf) 



i=3 



(38) 
» * 

-I{x s ;y d \h sr ,h sd ,h rd ) (39) 



We therefore wish to find 

C ^ C wors t = inf 

P»r(')>P» d (-).P z d(-) 

sup — 7(x s ;y d ,y^|/i sr , h s d,h r d)- (30) 
Px s (-) m 

The following result provides C wors t- 

Theorem 1: A lower bound on the capacity of AF scheme 
is given by 



where in d39l l we removed the dependence on i without loss of 
generality. Note that y is defined in ( f37T >. Now we can calculate 
the worst-case capacity by proving that Gaussian distribution 
for z r , z d , and z d provides the worst case. Techniques similar 
to that in [11] are employed. Any set of particular distributions 
for z r , z d , and z r d yields an upper bound on the worst case. Let 
us choose z r , z d , and z r d to be zero mean complex Gaussian 
distributed. Then as in [1], 



r 



< 



2m 



m — 2 
2m 



~E War E Wsd E Wri 



log (l+g(5 s ,P s ,a, 



-E log dct (I + (E(\x 8 \ 2 )AA r )(BE[zz^]B^)- 1 ) 

(40) 



f[g(8 s ,P s , <r sr ,w sr ), g{S r , P r , a rd , w rd )] 



sd, W sdj 



(31) 



where w sr ~ CA/"(0, 1), w sd ~ CN(0, 1), w rd ~ CAf(0, 1), 
and f(x.y) = xy/(l + x + y). Furthermore g(a,b,c,d) is 
defined as 



l{a,b,c, d) 



2a(l - a)ra 2 &V 



2(1 - a)mbc 2 N + (m - 2)(c 2 amb + N )N 



\d\\ 
(32) 



Proof : For better illustration, we rewrite the channel input- 
output relationships (1181 1. ( fT9l , and ( 120b for each symbol: 



y r [i] = h sr x s [i] + z r [i], 
yd[i] = h sd x s [i] + z d [i], 
for £ = 3,4,...,(m-2)/2 + 2, and 

V d [i] = KdX r \i] + z r d [i], 



(33) 
(34) 

(35) 



where the expectation is with respect to the fading estimates. 
To obtain a lower bound, we compute the mutual information 
for the channel d37l >. assuming that x s is a zero-mean complex 
Gaussian with variance i?(|a; s | 2 ), but the distributions of noise 
components z r , z d , and z r d are arbitrary. Thus, 

I{x s ;y d ; \h sr ,h sdl h rd ) = h(x s \h sr ,h sd ,h rd ) - h(x s \y d , h sr , h sdl h rd ) 

^ log7re_E(|2: s | 2 ) - log ire var(x s \y d , h sri h s d, h rd ). (41) 

From [11], we know that 

var(x s \y d ,h sri h sd ,h rd ) ^ E (x s - x s ){x s - x s )^ \h sr , h sd , h rd 

(42) 

for any estimate x s given y d , h sri h sd , and h rd . If we sub- 
stitute the LMMSE estimate x s — R xy Ry 1 fd into d4Tb and 
\, we obtain 



I(x s ;y d \h sr ,h sdl h rd )>Elogdet (I + (EWx^AA^BE^B*)- 1 ) 
As a result, we can easily see that 



for i = 3 + (m - 2)/2, (m-2) + 2. 

In AF, the signals received and transmitted by the relay have 

following relation: 



2m 



-£logdet (1+ {E[\x s \ 2 ]AA t ){BE[zz^}B'')- 1 ) 



(43) 



x r [i\=0y r [i-{m-2)/2], 



(36) 



From (l40l > and d43l . we have 

m-2 



£[|ar r 



2m 



.Blogdet (1+ {E[\x s \ 2 ]AA i )(BE[zz^]B J i 



\h sr \ 2 E[\x s \ 2 ]+E[\z r \ 2 ] 
Now, we can write the channel in the vector form 



Vd[i\ 

y d [i+(m-2)/2 



(44) 

Now combining d30l ), ((38) and ( 144b . and using the results ( T2"Tb . 
( 122b , d27b . ( |28T i and (|29j>, we obtain the following capacity 
lower bound 



1 

KdP 1 




_ m-2 

^worst - Eu u 



2 in 



log 1 



2(l-i,)mP, 
m-2 



isd | 



2(l-(5 s )mP 3 2 
m-2 °h. 



Nn 



(37) 



+/ 



2(l-S s )mP s 
m-2 



2(l-(5 r )mP r I £, . 2 



2(l-i,)mf, 2 
m-2 °)L 



N .> 20-Mm^ , 




(45) 



'Here we use the property that det(I + AB) = det(I + BA) 
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Intuitively, we may see the lower-bound as the case in which 
the estimation error is completely detrimental. By substituting 
(f6b-(fTTbinto d45l l and normalizing, we can rewrite the capacity 
as in (ED □ 

IV. A Capacity Lower-bound for DF 

In DF, there usually are two different coding approaches [2], 
namely repetition coding and parallel channel coding. We first 
consider repetition channel coding. For this case, an achievable 
rate is 

Irdf = — nim\l(x s -,y r \h sr ) 7 I(x s ;y d ,y r d \h s d,'h rd )\ . 
ml ) 

(46) 

Using this expression, we arrive to the following result. 

Theorem 2: An achievable rate expression for DF with 
repetition channel coding is given by 



where 

h = 



2m 



log 



(l + g(S s ,P s ,a s 



(47) 



(48) 



m-2 

2m Bd 



log + g(S s ,P s , a sdl w sd ) 
+ g(8r,Pr,o- r d,w rd ) 



(49) 



where g(.) is defined in d32l . 

Proof : As described in [11], we can obtain the worst-case 
mutual information for the first term in d46i i by proving that 
Gaussian distributed z r is the worst case. This gives us I\. 
In repetition coding, after successfully decoding the source 
information, the relay transmits the same codeword as the 
source. As a result, we can rewrite the data transmission with 
regard to the second mutual information in ( |46| > as 



y r d [i+(m- 



2)/2] 



Z 



h sd 

h rd (3 

A 

z d [i] 
dl i + (m- 



x s \i\ + 



2)/2] 



(50) 



z[,] 



In repetition coding 







E[\x s \ 2 Y 



(51) 



From ( l50l l. it is clear that the knowledge of h sr is not required 
at the destination. We can easily see that (T50b is a simpler 
expression than what we have in the AF case, therefore we can 
adopt the same methods as described in Section 3 to show that 
Gaussian noise is the worst case which gives li- The resulting 
lower bound capacity is 



st = mm 



{h,h}, 



where 

j _ m-2 
h 1m- E ^ 



log 1 + 



2(l-i>P, 
m-2 



2(l-i>f, 2 
m-2 a h s 



, (53) 



m-2 
h ~2m~ h 



log 1 



2(l-S s )mP s 
m-2 



\h sd \ 



2(l-i r )mf r 
m-2 



h>rd\ 



2(l-6 r )mP r 2 
m-2 a h r 



(54) 



Again by substituting (t6l>-dTTb into (|52b-(l54T> and normalizing, 
we obtain Theorem [2] □ 
If parallel channel coding is employed [2], then we have, 



iPDF 



— min {/(x s ;y r \h sr ), /(x s ; y d \ h sd ) + 7(x r ; y r d \ h rd ) X 
ml > 

(55) 



Again it can easily be shown that the worst case is experienced 
when z r ,z d , and z r d are Gaussian distributed. The resulting 
achievable rate is given in the following result. 

Theorem 3: An achievable rate expression for DF with 
parallel channel coding is 



where 



h = 



m-2 
2 m 



-E„ 



Iworst = min{/i,/ 2 }, 

log + g(S Sl P s ,a sr ,w sr )^j 



(56) 



(57) 



h = 



m-2 
2m 



~Ew sd Ew T d 



log (l + g(S s , P s ,a sd ,w sd )j 

+ log (l + g(6 r , P r , a rd , w rd )j 
where g(.) is defined in ( [32l . 

V. Optimal Resource Allocation 

We first study how much power should be allocated for 
channel training. In AF, it can be seen that 8 r appears 
only in g(S r , P r ,a rd ,w rd ) in the achievable rate expression 
PTT l. Since f(x,y) = . is a monotonically increasing 



(58) 



□ 



function of y for fixed x, ( BIT ) is maximized by maximizing 
g(S r ,P r ,a rd ,w rd ). We can maximize g(5 ri P r ,a rd ,w rd ) by 
maximizing the coefficient of the random variable |u; r( j| 2 , an d 
the optimal 5 r is given by the expression in (1591 . Optimizing 
5 S is more complicated as it is related to all the terms in 
(l3Tl l. and hence obtaining an analytical solution is unlikely. 
A suboptimal solution is to maximize g(5 s , P s ,cr sd ,w sd ) and 
g(S s , P s , cr S r, w sr ) seperately, and obtain two solutions ( 5* tt 1 f>opt 



and 5 s s U 2° pt , respectively. Note that expressions for J*" 1 °° :f ' c and 
^subopt ^ exact jy t jj e same as that in (|59i > with P r and a rd re- 
placed by P s , and a sd and a sr , respectively. When the source- 
relay channel is better than the source-destination channel, 

subopt 



subopt 



(52) g(S s , P s , cr sr , w sr ) is a more dominant factor and 5 S S U 2 ° v is 



4 



s° pt 



1 -AmP r a 2 rd - 2mN + 47V + 2^/-Am 2 P 2 a^ d - 2m 2 P r a 2 d N Q + m 2 N 2 - AmN 2 + 4JVg + 2m 3 P?(T$ d + m 3 P r a 2 d N Q 

2 -4mP r of , + m 2 P r a 2 

(59) 
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Fig. 1. <5 r vs. cr r( j for different values of P r when m = 50. 



good choice for training power allocation. Otherwise, 5 S S \ opt 
might be preferred. For DF, similar results and discussions 
apply. For instance, the optimal S r has the same expression 
as that in d59l >. Figure Q] plots the optimal S r as a function of 
a r d for different relay power constraints P r when m = 50. 
It is observed in all cases that the allocated training power 
decreases and convereges to a certain value with improving 
channel quality. 

In certain cases, source and relay are subject to a total 
power constraint. Here, we introduce the power allocation 
coefficient 9, and total power constraint P. P s and P r have 
the following relations: P s = OP, P r = (1 - 6)P, and 
P s + P r = P. Next, we investigate how different values of 
9, and hence different power allocation strategies, affect the 
achievable rates. An analytical result for 9 that maximizes the 
achievable rates is difficult to obtain. Therefore, we resort to 
numerical analysis. First, we consider the AF. The parameters 
we choose are P = 100, iVo = 1,5 S = 0.1, S r = 0.1. 
Fig0 plots the capacity lower bound (l3T1 l as a function of 
9 for different channel conditions, i.e., different values of 
&sr,&rd, and a s d- We observe that the best performance is 
achieved when 9 « 0.6 and a s d = 1, cr sr = 4, a r d = 4 which 
indicates that both source-relay and relay-destination channels 
are favorable. When a s d = l,c S r = 2,cr rc ; = 1, and hence 
the relay-destination and source-relay channels are not much 
better than the source-destination channel, the optimal value 
of 9 is close to 1 and there is only little to be gained with 
cooperation. 

Figs. |3]and|4]plot the DF achievable rates as a function of 9 
with the same parameters as in the AF case. Hence, the total 
power is P = 100. It is seen that paralel coding achieves a 
better performance compared to that of repetition coding. In 
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Fig. 2. AF achievable rate vs. 9. P = 100. (1) the dashed line cr a[ j = 
1, <T sr = 10, cr r d = 2; (2) the dotted line cr s[ j = 1, a sr = 6, <r r£ j = 3; 
(3) the dashdot line <r sd = l,cr sr = 4, <r r[ j = 4;(4) the solid line cr s( j = 
1, (T S r = 2, cr r[J = 1 



parallel coding DF, we observe that unless a sr is high and 
hence the source-relay channel is strong, the optimal value of 
9 is close to 1 and relay is allocated small power. 

We consider in this paper that there is a cost associated 
with cooperation. This cost is the power and time dedicated to 
learn relay-destination channel. This cost is more pronounced 
in the presence of a low total power constraint. Figs. [5] |6] 
and [7] plot the achievable rates when P = 1. We can see that 
DF have a better performance than AF at low power levels. 
Generally, cooperation gives more gains in the low power 
regime. However, as indicated by the solid-lined curves, if 
the quality of the source-relay and relay-destination channels 
is comparable to that of the source-destination channel, there 
is little or no gain through cooperation. 
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Fig. 5. AF achievable rate vs. 8. P = 1. 1) the dashed line cr sd = 1, <r sr = 
10, oy d = 2; (2) the dotted line cr sd = 1, cr sr = 6, oy d = 3; (3)the dashdot 
line u sd = 1, (Tar = 4, a rd = 4;(4) the solid line a ad = 1, Osr = 2, cr rd 
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